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EXECUTIVE  SUMMARY 


Vorticity  theory  has  been  successfully  used  to  theoretically  predict  and  explain 
lift  forces  on  wings,  ground  effects  during  aircraft  takeoff  and  landing,  and 
minimum  separation  standards  for  avoiding  trailing  vortices  from  large  aircraft. 
Although  aircraft  wings  are  relatively  large  areas  of  vorticity  generation, 
in  many  applications  these  regions  can  be  treated  as  vortex  filaments.  While 
aircraft  vorticity  generation  is  ultimately  caused  by  frictional  effects,  vortic¬ 
ity  in  fires  is  often  caused  by  pressure  and  buoyancy  forces.  If  a  fire  plume 
is  viewed  as  a  region  of  vorticity  production  somewhat  analogous  to  a  lifting 
surface,  the  possibility  of  approximating  flow  effects  by  a  vortex  filament  arises 
as  a  reasonable  approach. 

In  this  analysis,  two  classical  vortex  equations  (Bjerknes  theorem  and  the  Biot- 
Savart  law)  are  used  to  characterize  the  vorticity  produced  by  a  burning  wall. 

The  vorticity,  which  is  dispersed  throughout  the  flame,  is  idealized  by  a  single 
vortex  filament  located  in  the  upper  half  of  the  visible  flame.  By  employing  a 
mirror  vortex  element  behind  the  wall,  the  air  induced  into  the  flame  can  be 
calculated.  This  influx  of  air,  coupled  with  the  stoichiometric  requirements  for 
complete  combustion  of  the  fuel  evolved  from  the  wall,  leads  to  an  analytic 
relationship  for  the  size  of  the  flame  based  on  either  the  mass  loss  rate  of  the 
wall  or  the  energy  release  rate. 

The  final  relationship  agrees  well  with  theoretical  and  experimental  results 
recently  published.  The  use  of  vortices  to  describe  flow  fields  caused  by  fires 
offers  a  powerful  method  to  analyze  specific  fire  problems  that  might  be  extremely 
difficult  to  attack  by  other  techniques. 
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INTRODUCTION 


The  purpose  of  the  effort  is  the  development  of  a  simply  derived  relationship 
between  material  burning  rate  and  flame  height  on  walls. 


Background . 


Fire  growth  in  aircraft  cabins  involves  material  ignition,  burning  rates,  and  flame 
spread  rates  along  with  attendant  effects  on  the  gas  phase  environment  within  the 
enclosure.  The  fire  plumes  from  burning  materials  are  additionally  significant 
because  of  their  potential  for  directing  heat  at  nearby  or  adjacent  materials 
and  causing  fire  involvement  of  additional  fire  load.  For  instance,  if  the  visible 
fire  plume  of  a  burning  seat  impinged  on  the  ceiling,  flame  spread  on  the  ceiling 
surface  could  start,  provided  that  necessary  ignition  criteria  were  met.  Similarly, 
if  a  section  of  wall  were  burning,  the  height  of  the  fire  plume  would  dictate  how 
much  material  and  possibly  what  type  materials  were  exposed  to  direct  flame.  The 
height  of  wall  flames  has  direct  implications  for  upward  fire  spread  and  implicit 
bearing  on  material  selection  and  spacing.  These  implications  have  been  addressed 
through  analysis  and  experiment  (references  1,  2,  3,  and  4).  The  analyses  have 
been  done  with  dimensional  analysis  as  well  as  integral  methods  employing  boundary 
layer  type  equations.  These  approaches  involve  use  of  an  entrainment  coefficient 
to  account  for  the  introduction  of  enough  oxygen  to  support  combustion  of  the 
flammable  volatiles  from  the  pyrolysis  zone  of  the  burning  wall.  An  example  of  the 
type  relationship  developed  from  such  approaches  (reference  1)  is 
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where  Q  is  the  energy  release  rate,  Cp  is  the  specific  heat  of  air,  T0  and  PQ  are 
the  ambient  gas  temperature  and  density,  respectively,  g  is  the  gravitational  con¬ 
stant,  and  yf  is  the  flame  height.  Similar  equations  involve  a  different  coeffi¬ 
cient  from  4.2,  radiation  correction  factors,  as  well  as  different  arrays  of  terms 
used  to  factor  the  heat  release  rate.  The  striking  common  denominator  to  the 
experimental  data  and  theoretical  predictions  for  turbulent  wall  fires  is  that  the 
flame  height  is  proportional  to  the  two-thirds  power  of  the  heat  release  rate 
term. 


While  the  utility  of  entrainment  coefficients  is  well  established  for  turbulent 
jets,  wakes,  and  plumes;  there  is  another  approach  that  could  be  employed  in 
analysis  of  a  burning  wall.  The  alternate  approach  involves  viewing  the  flame  as 
a  region  of  vorticity  generation.  Once  the  vorticity  generation  is  modeled,  the 
air  required  for  combustion  will  be  provided  by  induction  rather  than  by  a  turbulent 
entrainment  mechanism. 


Objective. 


The  objective  of  this  analysis  is  the  use  of  classical  vorticity  theorems  to  estab¬ 
lish  the  relationship  between  flame  height  and  burning  rate  for  wall  fires. 
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Vorticity  is  a  defined  quantity  describing  rotational  movements  in  fluids.  While 
vorticity  is  most  often  associated  with  viscous  forces  in  boundary  layers,  it  can 
also  be  produced  by  pressure  and  density  gradients  in  fluids.  While  vorticity 
itself  is  defined  as  the  curl  of  the  velocity  vector,  an  additional  measure  of 
vorticity  is  the  fluid  dynamic  circulation  around  a  closed  loop  in  a  fluid.  This 
is  the  line  integral  of  the  vector  dot  product  of  velocity  and  the  loop  vector 
element.  The  starting  point  for  this  development  is  Bjerknes  Theorem  (reference 
5)  which  is  stated  as  follows: 


dr  -X  '(*)  x  (vp)  dA  (2) 

where  T  is  the  circulation,  and  the  right-hand  side  represents  a  surface  integral 
of  the  vector  cross  product  of  density  and  pressure  gradients.  In  a  wall  flame, 
the  pressure  gradient  is  vertical  and  the  density  gradient  is  horizontal.  Figure  1 
snows  a  schematic  of  an  idealized  wall  flame  and  the  corresponding  pressure  and 
density  gradients.  To  the  extent  that  these  gradients  are  perpendicular,  the 
vector  cross  product  will  give  a  unity  value  to  the  trigonometric  sine  term.  For 
the  wall  flame,  equation  2  can  be  restated  as: 


d 

dt 


dx  dy 


(3) 


where  J  is  the  width  of  the  integration  region  and  y  is  a  portion  of  the  flame 
height.  A  basic  assumption  now  made  is  that  the  flame  is  defined  by  some  average 
visible  colored  region  and  that  the  temperature  of  this  region  is  1800°  F  through¬ 
out.  In  this  way,  the  density  terms  in  equation  3  can  be  redefined  as  temperature 
terms  through  the  perfect  gas  law.  Figure  2  shows  an  integration  path  for  equation 
-J .  The  x-axis  integral  can  be  approximated  by  the  following  term  for  the  density 
gradient : 


JL  <tf  ~  ta) 

P5 


where  Tf  is  the  flame  temperature  and  TA  is  the  ambient  temperature  (equivalent 
to  T0  in  equation  1). 


Integrating  over  the  distance  6  the  term  comes  to 

R  (TF  ~  tA> 

P 

The  pressure  gradient  terra  is  that  caused  by  gravity  or  PA  g  which  can  be  written 
as  gp/RTA.  Thus,  the  integral  of  equation  3  on  the  right-hand  side  is 


If.  ~  iA 

TA 


g  y 


Assuming  a  steady  flame, 


(4) 


the  left-hand  side  of  equation  3  can  be  redefined  as: 


Thus 


V  dy  =  V2 


equation  3  simplifies 


to : 


v2  =  If  ~  -ta  s  y 

TA 


(5) 
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This  is  Decause  shrinking  6  to  zero  removes  the  line  integrals  in  the  x  -  direction, 
and  outside  the  flame  the  y  velocity  vector  is  approximated  as  zero.  Thu6,  Bjerknes 
theorem  leads  to  the  relationship  that  the  velocity  of  gases  in  the  flame  are  pro¬ 
portional  to  the  square  root  of  the  distance  from  the  base  of  the  flame.  Since  the 
circulation  around  a  loop  is  identical  to  the  total  vorticity  of  the  area  bounded 
by  the  loop,  equation  5  shows  that  the  vorticity  increases  higher  into  the  flame. 

At  this  point,  it  is  convenient  to  amass  this  dispersed  vorticity  into  a  single 
vortex  filament.  The  vertical  location  of  the  vortex  will  be  taken  at  the  point 
where  half  the  vorticity  of  the  flame  is  above  it.  Since  the  circulation  integral 
goes  like  the  integral  of  the  square  root  of  y,  the  location  of  the  vortex,  yv  will 
be  taken  when 

y  '  (6) 

v  2 

or  wnen  yf  is  1.59  yv.  Going  back  to  equation  5  and  integrating  the  velocity 
through  the  entire  flame,  the  total  vortex  strengh  at  yv  can  be  taken  as 

r  -  | IaTT  yf  3/2  (7) 

This  vortex  will  &e  located  at  a  close  but  unspecified  distance,  $v,  from  the 
burning  wall. 

At  this  point,  the  model  of  the  vortex  and  its  strength  does  not  provide  any  rela¬ 
tionship  to  the  burning  rate  of  the  wall.  To  develop  this  relation,  two  further 
relationships  are  needed.  First,  it  will  be  assumed  that  the  flame  tip  occurs  when 
enough  oxygen  has  been  brought  into  the  flame  to  combust  the  burning  wall  volatiles 
at  a  stoichiometric  ratio.  Second,  the  vortex  strength  must  be  used  to  find  the 
air  induced  into  the  flame.  The  first  of  these  can  be  stated  as: 


(?)6  4 


Tne  second  requires  a  measure  of  the  total  flow  induced  by  the  vortex. 

Figure  3  shows  the  type  magnetic  field  produced  between  two  buss  bars.  The  magnetic 
field  lines  are  the  analog  to  the  velocity  streamlines  in  vortex  theory.  From  the 
illustration  in  this  case,  all  inflow  ends  at  the  midsection  of  the  buss  bar. 

From  there,  outflow  begins.  In  this  vein,  the  gas  flow  between  the  vortex  location 
and  the  wall  will  be  determined  as  shown  in  figure  4.  That  is, 

"A*  V,  »v  <9> 

Figure  4  also  shows  a  mirror  vortex  on  the  other  side  of  the  wall.  This  vortex  is 
used  to  create  a  flow  field  preventing  velocity  components  normal  to  the  wall.  To 
find  the  flow  field  induced  by  these  vortices  of  opposite  sense,  the  Biot  Sarvart 
law  is  used  as  described  in  reference  6. 

V  «  T  /  ji_jc_r 

/  -3 


The  induced  velocity  from  each  filament  at  the  wall  will  consequently  be 

T  (11) 

V  ■  - 

1  2  w  <5 ,, 
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If  this  is  taken  a6  the  velocity  induced  in  the  flame  at  this  point,  then  the  volu¬ 
metric  flow  will  be  given  as 


V  0 
I  v 


'Employing  equation  7  and  12  and  multiplying  flame  density, 


PF  Vj  5 v 


=  OF  J.  \(Sl£  -  Ta)  g  yF 


Jsing  equation  8,  equation  13  becomes 


VJs  mf 


=  OF  J.  \  (Tf  -  TA)  g  yf 


Equation  14  shows  that  the  flame  height  is  proportional  to  the  burning  rate  (or 
energy  release  rate)  to  the  two-thirds  power  as  shown  in  equation  1. 

DISCUSSION 

Tne  form  taken  by  the  results  of  a  vorticity  analysis  differs  from  the  forms 
yielded  in  other  analyses,  primarily  from  the  fact  that  the  specific  heat  and  the 
entrainment  coefficient  do  not  appear  at  all  in  the  vorticity  analysis.  For 
comparison  with  equation  1,  which  has  the  parameters  used  for  data  analysis  in  this 
area,  some  reworking  of  equation  14  is  necessary. 

Rearranging  terms  results  in 


HL  3r 


"  1  Ta  ~~l 

\  ^TF“TA)  g 


Assuming  that  the  mass  loss  rate,  m j,  can  be  converted  to  energy  release  rate,  E, 
by  nui c i pi icat ion  by  some  heat  of  combustion,  AHc , 

E  -  n.;  A  He  (16) 


tne a  equation  13  can  be  written  as 


f  [\£  PMHc  2  \  Ta  (Tf_Ta)  g 


Inis  can  ho  converted  to  a  form  like  equation  1  by  multiplying  numerator  and  denomi¬ 
nator  Oy  the  rai->?ing  terms  to  get 


I  £  \  2^  \[~Tp  '  (Z  T  \l2/3r 

lTl  7aHcVt;Tt-^  (CpTMj  [■ 


Cp  TA  pa  /g-] 
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To  compare  equation  16  with  equation  1,  the  tirst  term  on  the  right  hand  side  has 
to  be  evaluated.  The  parameters  used  are  as  follows: 


AHc  =  10  k  J/g 


Cp  =  .24  cal/g  C 
Tf  =  1800°  F 

Ta  =  70°  F 


Using  these  numbers, 


E 

CP  Ta  P, 


_§ _ 1 

A  PA  fg  J 


( iy ; 


Given  the  uncertainties  in  flame  height  measurements  as  well  a3  the  different 
assumptions  used  in  a  vortex  approach,  the  close  agreement  between  equation  19  and 
equation  1  is  remarkable.  In  experiments  by  Hasemi  quoted  in  reference  3,  the 
coefficient  would  be  2.9  for  the  height  of  continuous  flame  and  6.1  for  the  upper 
most  flame  tips. 


The  common  finding  between  this  vorticity  approach  and  previous  studies  is  that 
E//ig  to  the  two-thirds  power  or  m//g^  to  the  two-thirds  power  is  the  most  important 
pair  of  terras  for  flame  height  correlations.  However,  the  development  of  the 
vortex  model  raises  questions  about  the  importance  of  other  factors  in  accurately 
modeling  flame  height.  The  vortex  model  induces  only  enough  air  to  combust  the 
pyrolyzing  material  while  the  entrainment  coefficient  used  in  reference  3  would 
involve  entraining  between  3  and  6.6  times  stoichiometric  air  into  the  flame.  This 
issue  could  be  resolved  by  experiment.  The  entrainment  theories  incorporate  the 
specific  heat  as  a  primary  parameter  while  the  vortex  theory  circumvents  this 
requirement.  Tne  reason  Cp  appears  in  the  right-hand  side  of  equation  19  is 
that  its  inverse  is  included  in  the  factor  3. 


Tne  vorticity  model  clearly  demonstrates  the  wall  flame  as  a  region  of  intense  vor¬ 
ticity  generation.  Although  this  vorticity  was  rolled  into  a  single  filament  for 
calcuLat  innal  purposes,  the  actual  flame  vorticity  is  of  necessity  dispersed  and 
disorganized.  In  actual  flames,  the  vorticity  cannot  roll  into  a  row  of  discrete 
vortices  because  the  wall  condition  would  place  mirror  vortices  directly  opposite 
to  them.  It  is  known  that  the  only  stable  situation  involves  staggered  vortices. 
Tms  vortex  instability  can  have  some  influence  on  wall  flame  tip  oscillation. 

Whue  the  vortex  model  dots  conveniently  explain  wall  flame  height  and  provide  in¬ 
sight  into  flame  stability,  the  major  weakness  of  the  model  is  that  it  does  not 
explicitly  give  a  flame  thickness.  Thus,  without  further  boundary  layer  type 
analyses,  the  vortex  model's  utility  is  limited  insofar  as  determining  heat 
transfer  to  the  wall  for  flame  spread  rate  prediction. 


CONCLL'S  iUWS 


uevelopment  of  a  vortex  mode f  to  predict  wall  flame  height  leads  to  the  following 
cone  1  as  ions : 

1.  Tt.e  wall  flame  height  is  proportional  to  either  the  mass  injection  rate  or 
tne  energy  release  rate  to  the  two -thirds  power  as  predicted  by  other  theories  and 
round  in  experiment. 

2.  The  coefficient  to  the  standard  heat  release  rate  array  of  terms  predicted 
by  the  vortex  model  is  well  within  experimental  uncertainty. 

.i .  An  imaginary  vortex  filament  can  satisfactorily  represent  the  dispersed 
vorticity  of  the  flame  for  calculations. 

4.  The  velocity  increases  with  square  root  of  height  as  the  flame  is  trav¬ 
ersed  vert ical lv . 

5.  The  vorticity  of  the  flame  causes  air  for  combustion  to  be  added  by  an 
induction  mechanism. 
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